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Abstract 



Domain structures appear in ultratliin magnetic films, if an easy-axis anisotropy, 

favouring the direction perpendicular to ttie film, overcomes the shape anisotropy. We 

1^ \ will study the nature of the magnetic ground state and how the different parameters 

r — ' influence the domain structure. Both Ising and Heisenberg spins will be considered: in 

^^ ■ the latter case, a reorientation phase transition may take place. Domains are expected 

^\i i to appear just before such transition. Domain structures at finite temperature will 

be described in terms of a two dimensional array of domain walls, whose behaviour 

resembles a liquid crystal. 



g '. 1 Introduction 

pH ■ A three dimensional ferromagnetic sample, below its Curie temperature, generally has no 

O . net magnetization because —as suggested by Weiss [H and discovered by Bitter |^— it is 

made up of contiguous regions where the magnetization is oriented along different directions. 

These regions (typical size ~ 1 ^ 10/xm) are called magnetic domains, and are separated 

r> I by domain walls (typical size ~ 10^ -i- lO^A), where the magnetization rotates continuously 

C^ ' from one direction to another. 

The magnetization inside each domain is due to the ferromagnetic coupling (Tiex) be- 
tween neighboring spins, which aligns them parallely. The existence of domains is due to the 
dipolar energy (Tidip) : the true "magnetic" (and long-range) interaction between spins. This 
energy is also called "magnetostatic energy" or "shape anisotropy": in fact, its minimiza- 
tion strongly depends on the shape of the sample. Finally, domains are not magnetized in 
random directions, but they prefer certain directions, depending on the crystal-axes, whose 
existence is due to the so-called magnetocrystalline anisotropies (Tian) • The actual magnetic 
structure of a system depends on the minimization of the total energy Ti. = Tiox+'^dip+'^an- 
In this paper, we will be interested in ultrathin magnetic films: by "ultrathin" we mean 
a thickness smaller than, or of order of A^ = 10 atomic monolayers. Such small thicknesses 
can be currently obtained through growth by Molecular Beam Epitaxy 0. Our aim is to 
give a microscopic description of domain structures in such two-dimensional systems, both 
at zero temperature (Section ||) and at finite temperature (Section H). Section is devoted 
to the dependence of the domain size L on the film thickness A^. A comparison with Kittel's 
theory —valid for much thicker films (A^ > L)— will be done. 



Ultrathin films can undergo a reorientation transition, with the magnetization passing 
from being perpendicular to the film to being parallel to it. Since the condition for the 
occurence of this transition coincides with the condition for having domains not unphysically 
large, the interplay between reorientation transition and magnetic domains will be treated 
in Section |^. 

2 Domain structures at zero temperature 

The aim of the present section will be the study of the classical ground state of a ferromag- 
netic ultrathin film, in presence of dipolar interactions and of a uniaxial anisotropy. Here 
and in the following, both the spin modulus \S\ and the lattice constant a of the cubic 
crystal will be put equal to one. 

The hamiltonian of a single monolayer will generally be assumed to have the form 

H = Hex + Han + ^dip (1) 

where the first term 

n,x = ^Jdx{vsy (2) 

is the Heisenberg interaction, in the continuum approximation. The second term 

n.n = -^ j dxSl{x) (3) 

is a single-ion anisotropy caused by the crystal field; it appears because, in a thin film, 
the z direction (perpendicular to the film) is not equivalent to the in-plane directions. 
The constant K may be negative (easy-plane effect), but in this case no domain structure 
appears. Therefore, only the opposite case {K > 0) will be considered. Finally, the third 
term 

^dip = n^Z + ^iip (4) 

is the dipolar interaction, which can conveniently be split into a rotationally invariant part 

<=l^/#^S(f).S(f') (6) 

and an anisotropic part 



dip o / I -* -•/ n * 



which couples the spin space to the "real" space. In addition to the constant fi, we have 
introduced the notation: z7 = (x — x')/\x — x'\. 

Since a spin doesn't interact with itself, the integration must exclude a region around 
X = x' . This problem doesn't exist if a discrete calculation replaces the integration (see 
below), or in the calculations relevant for domain structures, where x and x' are always far 
away (because belonging to different domains). 

The term TH/^ represents a long-range, antiferromagnetic coupling, while 'h}^ introduces 
an easy-plane anisotropy. In fact, this term vanishes for a perpendicular coUinear state 



(hereafter called _L state), because S ■ u = 0, whilst for a parallel coUinear state (hereafter 
called II state), it has the negative value: 

8tt^ |x|3 ^' 

per spin. So, the || configuration is energetically favoured with respect to the X one; 
furthermore, since only the in-plane component S\\ of the spin contributes to it, for a generic 
coUinear configuration forming an angle 6 with the z axis, {S ■ u)"^ rewrites as (5*11 ■ u)"^ 
and therefore expression (|^) simply needs to be multiplied by sin^ 9. So, for a coUinear 
configuration, TCdip is found to reduce —up to a constant— to the anisotropic expression: 

n'^^^ = - ^ sin^ e f dx=^ f dxSl (8) 

Sometimes, instead of Vt we will use the quantity g = Q/lGn, which appears in the 
discrete version of dipolar interaction: 

"Hdip = 9^^— Si- Sj - 3{Si ■ u){Sj ■ v) 

Because of the competition between T^^ip (^ weak long-range antiferromagnetic (AFM) 
coupling) and ?^cx (a strong short-range ferromagnetic (FM) coupling), a domain structure 
arises, for any value of the dipolar coupling; furthermore, the competition between Tidip 
(an effective easy-plane anisotropy) and Tian (an easy-axis anisotropy) can produce a re- 
orientation phase transition (RPT) of the magnetization, from the _L state to the || state. 
To prove the former statement, it is sufficient to show that the _L state is destabilized by 
the introduction of an Ising domain wall. So, if we invert in the negative z direction all 
the spins in the half-plane x > 0, we have to compare the cost of the domain wall (2J 
per unit length) with the energy gain, due to dipolar interaction. Such gain is obtained by 
summing up {—2g/\x — x'\^) on all the spins (x') of the half-plane x > and the spins (x) 
of the half-plane x < 0, with y = constant (because we need the gain per unit length). The 
first sum on x' corresponds to two spatial integrations, giving a quantity of order — 25f/|x|. 
Therefore, the third sum diverges logarithmically at large distances^ whatever the value of 
g, thus proving that the dipolar gain always prevails on the exchange loss. 

The previous rough calculation also allows to have a semi-quantitative expression for 
the domain size of a stripe domain structure. In fact, the third integration on x will be 
limited up to the domain size L, by giving an energy gain of order —glnL which must be 
approximately equal to the domain wall energy (per unit length): glnL ^ Edw So: 

L ^ exp{Ednj/g) . (9) 

For an Ising system, E^w = 2J, so that (at zero temperature: T = 0) the domain size is 
exponentially large, unless J and g are of the same order. In Section ^ we will see how the 
previous result is modified, if thermal effects are taken into account. 

For an Heisenberg system, the domain wall energy has the form [Q: E^w = 2vAJ, where 
A is the total effective easy-axis anisotropy (see Eq. (^)): X = K — Q. Therefore, domains 
of "small" extension can be obtained even ii J ^ K, Q, provided that /C ~ fi. 

-'^At small distances, the lattice constant is the natural cut-off. 



Our calculations will concern a square lattice. A triangular lattice gives rise to frus- 
tration, in presence of an antiferromagnetic exchange coupling. Indeed, dipolar coupling is 
an antiferromagnetic interaction, in the Ising limit; anyway, it is a long-range one: so, the 
induced frustration exists even in a square lattice. We expect that the picture discussed in 
the rest of this Section is not qualitatively modified for a triangular lattice, if L ^ 1. 

2.1 Domain walls 

The size of a domain wall is given by the balance between the effective easy-axis anisotropy, 
and the exchange energy. If the long-range character of Tidip is neglected and only its easy- 
plane effect (D is retained, it is possible to determine analytically the magnetization profile 
inside the wall, and the domain wall extension. 

If the wall is parallel to the y axis (as it will be always supposed in the following), in 
the continuum approximation the domain wall energy (per unit length) is 



'Hdw = - I dx 



J{VSf-\Sl , (10) 



which must be minimized with respect to S{x). 

By using polar coordinates {S = {sm6cos(j),sm6sm(j),cosd)) we obtain 

n,^ = jY{A{dJ? + {d.<P?sm'e]-Xcos'e} . (ii) 

The 0-dependent term is minimized when is constant, e.g. = 7r/2 (Bloch wall) or 
= (Neel wall). On the other side, the minimization with respect to 6 gives 2d'^^6 = 

(A/J) sin(2^), whose solution gives 9{x) = 2arctan fe^/"^), where w = Jj/\- The corre- 
sponding expression of Edw is obtained from Eq. ([T0|) : Edw = J/ui + wX = 2\fXJ . 

The degeneracy among different types of wall is removed if dipolar interaction is correctly 
taken into account. It is immediately seen that the relevant interaction is the self-energy of 
the domain wall: in fact, let us consider the dipolar interaction between a spin 5*1 belonging 
to the wall, and a spin 5*2 belonging to a nieghbouring domain. Since 5*2 ■ i^ = {v = r/r, 
f^ being the vector joining the two spins) only the isotropic interaction g{Si ■ S2)/r^ is not 
vanishing. Anyway, the quantity 5*1 ■ 5*2 = cos 6i does not discriminate among different types 
of wall, because it only depends on the axial angle di of the spin 5*1, whilst walls differ in 
the polar angle 0i. 

The dipolar contribution to the domain wall energy therefore writes: 

J,d^p^9_^^ f. f S{x).S{x') {S{x)-r)iS{x')-r) \ 

where we have used the translational invariance in the y direction, and the following nota- 
tions are understood: r = {x — x',y) and 5'(x) = (sin6'(a;) cos0, sin6'(a;) sin0, cos6'(x)). The 
function 6{x) is any increasing monotonic function, such that ^(0) = and 6{w) = it. Our 
aim is to minimize Ej^ with respect to the polar angle 0. 
It is straightforward to obtain the following relations: 

s{x)-s{x') = cos{e{x)-e{x')) 

{S{x) ■ r) {S{x') ■ r) = sin 9{x) sin 9{x') [{x — x'Y cos^ + 1/^ sin^ (p + 2y{x — x') sin cos 0] 



By retaining only the 0-dependent terms, we have: 
^^.p_ 3^AA ^ sin^(a;)sin^fx' 
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The term proportional to sin 20 gives a zero contribution, because it is an odd function 
of y. So, we are left with the expression: 
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Generally speaking, the shape anisotropy due to dipolar interaction implies that two 
spins S{x), S{x') will prefer to align parallely to (x — x'). Therefore the sum T{w) simply 
reflects the shape anisotropy of a stripe: y'^ > {x ~ x'Y, where the bar means some average 
on the stripe. In other words, the dipolar energy of the domain wall is minimized if the 
spins point on the average in the (elongate) y direction. The positive factor sin^(a;) sin^(x') 
does not play a major role: it is present because the in-plane component of the spin varies 
inside the wall. In conclusion, T{w) is positive and Ej^ is minimal for = 7r/2, i.e. for a 
Bloch-type wall. 



2.2 The Ising model 

Once we have shown that the ferromagnetic configuration is destabilized by dipolar inter- 
actions, whatever the value of g, in this section we will study more closely the Ising model, 
in the entire range of values for the ratio g/ J . Let us start with the physically relevant case 
g/ J <^ 1, which corresponds to very large domains. 

By extending a little more the naive calculation made in the introduction, it is possible to 
determine the leading term in the expression of the dipolar energy of a stripe configuration. 
In order to use the following results also for the Heisenberg model, and for a larger number 
of atomic planes, we will suppose the two neighbouring stripes of size L to be separated by 
a domain wall of size w <^ L and to lie on different planes at distance c. So, the relevant 
integral is the following: 



/(c) = - 

whose general solution is: 
'4L - 2w 



rL—w rL r+co 

/ dx' dx dy 

Jo Jo J-oo 



Kc) 



c 



/2L-w\ , 
arctan ( I — In 



[{x + x'y + y^ + c^f^ 



(13) 
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In the case c = (corresponding to a single monolayer) or for w ^ c, we find /(c) ^ 
21n(L/2u'). In these cases, the dipolar energy gain with respect to the _L state is, for a 
couple of planes at distance c: 



7-i2 planes 
^dip 



.2<,/(c) = -4<,l„(|-) 



(15) 



For an Ising monolayer, c = and w = 1, the total energy per spin of a stripe (S) configu- 
ration with respect to the _L state is: 

E^(L) = ^-|ln(L/2). (16) 

However, the factor 2 inside the logarithm (which modifies, in a certain sense, the value of 
J: J ^ J + 2gln2) is not meaningful, because it contributes to the non leading term of the 
dipolar energy. Conversely, it is noteworthy that the leading term (proportional to g In L) 
is the one obtained by considering the interaction (per unit length) of a stripe of size L 
with a neighboring stripe of infinite size! Analogously, the same leading term appears if the 
neighbouring stripe is indeed a square of size L. This result therefore provides —in the same 
approximation— the expression E^[L) for the energy of a square network (C=chess-board) 
of size L: 

i?^(L) = ^-|ln(L)^2E^(L). (17) 

Therefore, not only the exchange cost of a square network is the double of a stripe structure 
with the same period (as the total length of domain walls is the double), but also the 
dipolar gain is the double. In very poor terms, the reason is that the infinite size of a stripe 
contributes only to the next-to-leading order term. Since in a square network each domain 
has four neighboring domains with opposite magnetization, whilst in a stripe configuration 
only two, we have the result E'^{L) = 2E^{L) (at the leading order in the dipolar energy!). 
The previous results have been found by Czech and Villain ^ (see their Eq. (12)), who 
have drawn the erroneous conclusion that the chess-board configuration has a lower energy 
than the stripe configuration. Indeed, the next-to-leading order terms modify the previous 
picture, as explained by Kaplan and Gehring @. Through a more careful evaluation of the 
dipolar sums, they find that E^[L) and E^{L) are modified as follows: 

^t?/ ^N 2 J — acq 4:0 , ^ 
E^{L) = — ^-^InL 

Ij Lj 



where as — 4.2 and ac — 0.3. By minimization with respect to L, we find 

Ls = exp 
Lc = exp 



2g 4 



^ + 1-^ 
2g 8 






From the previous expressions, it is apparent why neglecting terms of order [g/L) in 
the energy leads to mistakes. Furthermore, E^ < E^ if Lc/ Ls > 2, that is to say {as/4: — 
ac/8) > In 2, a relation which is satisfied. 

Now, we will consider the opposite limit of a large ratio g/J: dipolar interaction prevails 
over the FM exchange. In the following, the notations C-L and S-L will mean, respectively 
a square network of side L and a stripe configuration of period L. Since dipolar interaction 
is a (long-range) AFM coupling, we are led to conclude that the AFM state (i.e. the C-1 
state) is the lowest energy one: an explicit numerical calculation |0 shows that the C-1 
state is indeed the ground state. 



Such a state is maintained for a sufficiently weak J. So, a square network is preferred 
at (very) high g/ J , but stripe configurations are preferred at low g/ J . We can ask when 
the transition between square configurations and stripe ones takes place. According to 
numerical results 0, for J ~ 0.4(yf, the C-1 state is replaced by the S-1 state, and for 
increasing J we simply have an augmentation of the period L of the stripe configuration. 
Since analytical theories are not applicable here, let us give some simple arguments to justify 
why C-1 is replaced by S-1, rather than C-2. Let us compare the energies of the two possible 
configurations: 

E\\) = Ei^{l) + 2J. 

So, it is sufficient to compare their dipolar energies. Both the configurations have a vanishing 
nearest bond contribution to the dipolar energy (two spins up and two spins down), but in 
the S-1 state all the four nn-bonds are AFM, whilst in C-2 are half and half. For the 3^^^ 
nearest bonds, the picture is reversed, but the coupling is weaker, and up to the 6*^ nearest 
bonds the contributions are the same. This simple evaluation gives the correct result. 

We conclude this section by summing up the main results. Dipolar interactions favour 
an AFM coupling between spins; if we consider chess-board (C-L) and stripes (S-L) config- 
urations, we have that the gains of dipolar energy (with respect to the _L state) \E^^p{L)\ 
and |£'^p(L)| are decreasing functions of L and |i?^p(L)| > \E^^p{L)\. Nevertheless, when 
a FM exchange interaction is taken into account, the stripe configuration is energetically 
favoured, except for J/g < 0.4, when the ground state is the AFM one. 

2.3 The Heisenberg model 

— * 

In the previous section we have supposed that spins are "Ising spins" : S is always parallel 
to the z axis. Now spins are allowed to rotate, and the main difference is the energy of the 
domain wall, which is lowered with respect to the Ising case. This allows smallest domain 
sizes, even if —at T = 0— they will be seen to remain of macroscopic extent, except in a 
very narrow region of the parameter space. A second important feature of the ferromagnetic 
Heisenberg model, is that its excitations are spin-waves, which allow a local linear stability 
analysis. 

Some insight in the model is given by the study of the collinear configurations, the 
generic one forming an angle 6 with the z axis. Since the total effective uniaxial anisotropy 
is A = {K — Q), the _L state {6 = 0) will be the lowest energy one for K > Q. For K < Q 
the magnetization will lie in the plane [6 = tt/2), and for K = Q all the collinear states will 
be degenerate. 

The previous section has been devoted to show (for an Ising model) how the long-range 
character of the dipolar interaction "destabilizes" the _L state. This result is certainly 
not modified for an Heisenberg model, which simply allows for a much larger phase space. 
Anyway, there is no reason to think that an infinitesimal variation in the _L state will lead 
to a lower energy: in other words, the ± state may be locally stable and spin-waves do not 
destabilize it, i.e. it may be a metastable state. In fact, the dispersion curve writes: 

uj±{q) = {K-n) + Jq'^ = X + Jq\ (18) 

Dipolar interaction simply modifies the gap, because of its "easy-plane" effect. 



If now we turn to the || state {K < Q), a. first question is: Do dipolar interactions 
destabilize the || state? The answer is no: this case, domains are not energetically favoured. 
In the _L state, domains rise because of the competition between the FM Hex and the AFM 
Ti^jp. Anyway, in the || state, the anisotropic term Ti^ip (which vanishes almost everywhere 
in the _L state) changes the character of T^dip- Let us compare the in-plane FM and AFM 
configurations. 

The spin in the origin is choosen along the y axis: this way, all the other spins can be 
arranged in circular rings of increasing radium where spins are mutually parallel, and —in 
the FM state— always oriented along the y direction, or —in the AFM state— alternately 
oriented along y and —y. An in-plane domain structure may arise only if Tidip favours the 
AFM configuration with respect to the FM one. To check it, let us consider the dipolar 
interaction between the spin i in the origin and all the spins j on a given circular ring at 
distance r. The resulting sum is proportional to ±[1 — 3{{S ■ i^ij)^)j]/r^, where the sign -|- 
(resp. — ) holds for the FM (resp. AFM) configuration. Since the average (. . .)j on the spins 
j gives (1/2), the negative anisotropic part prevails on the first one, and the energy is lower 
for the FM state! Consequently, domains are not expected if the magnetization lies in the 
film plane. 

Anyway, the linear stability analysis of the || state reveals important features. In fact, 
it was shown by Maleev [^] in the seventies that dipolar interactions are able to stabilize 
the II state at finite temperature, because of their long-range character. This means that 
the dispersion curve is modified well beyond the introduction of some effective easy-plane 
anisotropy, which is unable to assure a finite magnetization. The correct expression for the 
dispersion curve is g, ^, |10[ {9^ is the in-plane angle between g and the magnetization): 



^fliq) 



{n-K)-^q + Jq' 



X 



—qcos %+ Jq 



(19) 



We have assigned a different meaning to the two terms, because now the Goldstone mode 
{uj{q = 0) = 0) exists for Uin, but not for ujout- Anyway, an out of plane excitation with 
a finite wavevector may go to zero and destabilize the || configuration. This happens for 
q = Q/AJ, when 1 > (K/Q) > 1 — {Q/16J). So, for a weak dipolar interaction such interval 
is extremely narrow. 

We can summarize the results of the spin-wave stability analysis of the FM configura- 
tions, as follows. The ± state is locally stable, and therefore it is metastable, because the 
insertion of a domain wall lowers the energy of the system. Domains are not expected in 
the in-plane state; nonetheless, the long-range character of Tidip destabilizes the || configu- 
ration, if Q is just a bit larger than K. The main result of the next section will be to link 
quantitatively the results of the previous analysis, with energetic calculations on domain 
structures. 

2.4 Yafet and Gyorgy's theory 

To my knowledge, the only microscopic study of domain strutures in a Heisenberg system, 
is due to Yafet and Gyorgy ||TT|. They study a stripe domain structure, characterized by 



two parameters: the domain size L and the domain wall size w. The wall is of Bloch-type, 
and the component 5*2 of the magnetization, along the domain wall, has a sinusoidal profile. 
As suggested by the same authors, we have added a third parameter which reveals essential 
to match the results of the spin-wave stability analysis: the maximal component s of the 
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Figure 1: (a) The maximal value s of the component \Sz\, inside each domain, and the 
domain size L (in logarithmic scale), (b) The energy gain of the domain structure, relative 
to the coUinear state (the _L one for / > 1 and the || one for / < 1). The smallest value of 
/ shown in (b) is fm, to which corresponds —in (a)— the vanishing of s. For both figures, 
2Ld = 25. Similar dependence for L{f) and (b) had already been given in Ref. ||TT |. 



magnetization in the z direction, inside each domain, s = means the || configuration, while 
s = 1 is a "standard" domain structure. The energy of a generic stripe domain structure is 
therefore a function of the three variational parameters {L,w,s). Its minimization, and a 
comparison with the energy of a generic coUinear state, will allow to determine the lowest 
energy configuration. 

In Fig. |I| we plot (a) L and s, as a function of the ratio / = K/Q, and (b) the energy gain 
of the domain structure. Two features are manifest: (i) It exists a minimal value /^ < 1, 
beyond which a domain structure is energetically favoured. As expected, the energy gain 
persists till / = oo and it is maximal for / = 1. (ii) When f = fm, spins lie in the plane 
(s = 0), and L = Lm = 4:7iLd, where Ld = J/Q is the dipolar length. Upon increasing /, 
the maximal component of Sz immediately saturates to s = 1, and the domain size L grows 
exponentially: unless dipolar coupling and uniaxial anisotropy nearly compensates, domains 
are unphysically large. How close to one / must be, in order to observe domains, it depends 
on the dipolar length. It is not possible to give here the details of the calculation, which 
follow Ref. [^ . Anyway, the introduction of the parameter s modifies the values of fm and 
Lm found by Yafet and Gyorgy, which now match the corresponding values, determined via 
the linear stability analysis of the || configuration |T^, I3\: fm = ^~ (fi/16J) and L^ = 7r/g 
(see below Eq. {^)). 

In conclusion, the following picture emerges: If K is fairly smaller than Q, the || config- 
uration is the ground state; when K is nearly equal to Q {K = Q — Q'^/16J) a soft mode 



appears, with a finite wavevector q = ir/Lm- In the new configuration, spins acquire a 
component S^ which is alternately positive and negative. Upon further increase of K, a 
domain structure clearly appears, with regions where 5*^ = ±1 are separated by domain 
walls of size w. These configurations are the true ground state, but for K > Q the _L state is 
locally stable. The domain size L increases exponentially and rapidly becomes larger than 
the sample size. 

The main difference between Ising and Heisenberg models is the form, and therefore the 
cost, of a domain wall: in the former case, the wall size is one lattice constant, and its energy 
per unit length is simply E^^"'^ = 2J. In the latter case, the wall can extend over many 
lattice constants and its cost is given by the balance between easy-axis anisotropy and 
exchange: £'^J**<=" ^''^ = 2J J{K — Vl). This expression correctly shows that the domain 



wall energy goes to zero, when K and Vt nearly compensate. Nevertheless, it has the 
drawback not to match E^^^^ when K -^ oo: this is because it derives from a continuum 
approximation, which fails to reproduce the discontinuity of an Ising domain wall. However, 
when the Ising limit is pertinent, domains are very large and probably unphysical. So, this 
limit is not really relevant (at least at T = 0). 

An approximate analytical expression for L can be found by noticing that the dipolar 
interaction between domains can be evaluated —in the Heisenberg model— by "removing" 
the domain wall fl^: so, we have Ising domains left, at distance w. By using Eq. ([l^ ) and 



the espression given in Section 2.1 for the energy of a domain wall, we have the followig 



formula for the energy per spin of a stripe domain configuration: 

^ Aw J 4gf , / L \ , , 

which must be minimized not only with respect to L, but also w. The calculation gives: 



w 



2[l + ^l+AH^La{f-l] 
L = 2wexp{J/2gw) . (21) 

This expression for L is valid also for the Ising model, where w = 1 The form of w proves 
the existence of a minimum value fm, for which the square root vanishes and L c::^ L^. 
The value of fm found via the approximate expressions given above is not the correct one: 
/m = 1 — l/16Ld, because Eq. (pO|) neglects the parameter s. 

3 Domain size and film thickness 

3.1 Kittel's theory and modifications 

The theories developped in the previous section apply to a single monolayer. Anyway, before 
extending them to a (thin) film of A^ planes, let us report Kittel's old results JlSl, which 
are valid for a thick film, whose thickness is larger than the domain size (A^ > L), and a 



recent modification |T6|, introduced to take into account the possible surface character of 
the anisotropy K. 

In Kittel's theory, the total energy [Fdom) of the magnetic slab is the sum of three con- 
tributions: the magnetostatic one (Fdip), the domain-wall energy (F^^), and the anisotropy 
term (Fan)- He considered both the possibility of fiux closure (Fig. ^a), and of an open 
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Figure 2: Domain structures in a thin film: fiux closure (a) and open fiux circuit (b). 

fiux circuit (Fig. ^), because he was interested in determing the lower energy configura- 
tion. However, the dependence L{N) is the same for the two cases. If closure domains are 
present, the contribution of the magnetostatic energy is neglected and the anisotropy term 
is proportional to the volume of the domains where the magnetization lies in the film plane. 
Finally, F^w depends on the extension of the domain walls. If a^i and cr^2 are the wall 
energies per unit area, respectively of the 90° and 180° walls, the total energy -FJ"^ (per 
unit of surface) of the closure domain structure, writes^: 



dom 



pi'^) I p{a) ^ 
dw 'an 



2a^i , fN 

\- cru>2 T tan a 

cos a \L 



+ K— tana 
2 



(22) 



Kittel considered the limit N ^ L and assumed a^i = o"„,2 and a = 45° 



The minimization of F!jZ„ gives L 



2cr^2N/Kt8ina and FJ"^ 



dom &^^^^ -^ — Y^""'2^v/jv ^aiiLx aiiu. ±^^^ — (2o-^i/ COS «) - 

aw2 tan a + yflKtaiuAxOwiN ■ So, the domain size is proportional to the square root of the 
film thickness. 



Such a result is found also in absence of closure domains: in such case F, 



(.b) 



dom 



.(ft) 
dw 



F^ ' 4- F 



(b) 
dipi 



because FJ^ = 0. Nevertheless, the relevant functional dependences are unchanged: Fj^ — 
<yw2N/ L and F^^^ ^ M^L {Ms being the magnetization of saturation), i.e. it has the same 
L-dependence of Fj^) . Consequently, the result L r- 



N still holds. 



Recently, the above calculation has been modified |T^ to take into account the possibility 
of a surface anisotropy, in an ultrathin film. In this limit, N <^ L, closure domains are 
supposed to make a small angle a with respect to the surface of the film, and the total 
anisotropy K is written as the sum of a volume {K^) and a surface (Kg) part: K = 
Ky + 2Ks/N. If A's(~ 1 erg/cm^) prevails over K^{^ lO^erg/cm^), the result by Kittel is 
modified as follow^: 

(23) 



2a.„oN 



Cr«,2 



-N . 



V Ktana y i^^tana 
So, for a surface anisotropy, Kittel's theory gives a linear increase of the domain size with 



the thickness. This result agreed with some experimental findings |17| on thin epitaxial 
Co/Au(lll) films. Nevertheless, the calculation in Ref. |TB| does not seem to be consistent. 



^The superscripts ^"-' and '^''^ refer, respectively to Fig. 2a and Fig. 2b. 
^In Ref. |lq|, the trigonometrical factor tana is missing. 
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In fact, it is noteworthy that below a critical thickness, Kittel's theory predicts a single- 
domain in-plane state. This result is obtained by comparing the energy Fy = 2Ks + K^N 
of a single-domain in-plane state with the energy of the closure domain state: 

Pdom = — ^ + '^'JKs(Jw2 tan a - cr^2 tan a . (24) 

cos a ^ 

If we write a^i ~ c^^„2 — cr«, and we neglect the last two terms with respect to the first one 

(because tana ^ 1), we find that the in-plane magnetized state is energetically favoured 

ifQ A^ < 2(0"^; — Ks)/ K^: a condition which is generally fulfilled in ultrathin films. 

A more rigorous consideration is the following. Expression (^2|) for the energy F^l^ of 

a closure domain structure should be minimized not only with respect to L, but also to the 

angle a. If daFj^ is calculated, we obtain 

daFt± = -^ [{KL/2 - 0-^2) - 2(7^1 sin a] = . (25) 

cos^ a 

However, it is easily seen that its solution: sin a = (i^L/2 — (j^2)/2o"^i does not correspond 
to a local minimum: the Jacobi matrix of the second derivatives has a positive and a negative 
eigenvalue, because F^'^ has a minimum with respect to L, but a maximum by varying 
a\ Indeed, this problem exists also for the original Kittel's theory, where a is arbitrarily 
put equal to 45°. A correct approach should probably add an a-dependent magnetostatic 
contribution F^^^. 

3.2 Ultrathin film limit: microscopic approach 

In the following, we will be concerned with the extension of the Ising and Heisenberg 
models to the case A^ > 1. Whilst the extension to A^ > 1 of the domain wall energy is 
straightforward for both models, the contribution of the dipolar interaction energy between 
different planes is not trivial, and some further approximation is necessary, in order to 
obtain tractable formulas. 

Let us start with the Heisenberg model, which is even simpler than the Ising one. The 



domain wall size w varies in the interval J J/K < w < J/g. So, the hypothesis w ^ N 
is reasonable, and in this limit the integral I{c) in Eq. ( |14|) does not depend on c. Such 
integral appears A^ times for the value c = and A^(A^ — l)/2 times for c 7^ 0, that is to say 
a total of A^(A^ + l)/2 times. Therefore, the expression (|20| ) for the energy of the domain 
structure is modified as follows: 

E,om{N) ^^ Xw ^ J7(A0 _ 2(Ar-H)^ ^^ /^x _ ^^^^ 



N L Lw L \2w, 

The factor 7(A^), whose importance will be elucidated very shortly, has been introduced 
to take into account a possible "long-range" character of the exchange interaction. If the 
range r^x of ?^ex is much smaller than A^, then 7(A^) = 1, but 7(A^) = 7*(A^) = (A^ -|- l)/2 
if Tex > N. Finally, the anisotropy term Han has been assumed to be of volume type. 

The A^-dependence of Edom can be reduced to the following transformations: J -^ 'y{N)J 
and g -^ '~f*{N)g. The domain wall size w has no trivial transformation, but in the limit 

L > w, w ^ ^J-f{N)w (see Eq. ^. So: 

I V-/{N) J 

L(A^) ^ •\/7(A^)u7ie^*('^) 29-1 , (27) 



*A11 the energies must be expressed in units where the lattice constant a — 1. 
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where wi refers to the monolayer (A^ = 1). Since J'y{N)/'y*{N) is always a decreasing 
function of N, a strong shrinkage of the domains is expected, even if Tiex is "long-range". 

Now, let us consider the Ising model. When we deal with the dipolar term, it is no more 
possible to suppose w ^ N, because now w = 1. The integral (|l3|) writes: 



/(c) ~ 2 + (2/c) arctan(l/c) + 2 ln(L/2Vl + c^) ~ 2 ln(L/2Vl + c^) , (28) 

if L ^ c. The total dipolar energy gain E^ip is obtained by summing up /(c) over c, for all 
the couples of planes. Kaplan and Gehring ||^ have used a continuum theory by resorting to 
the magnetic potential: this corresponds to the limit A^ ^ 1. In this case /(c) ~ 21n(L/2c) 
and 



N L L \N, 

The factor N inside the logarithm does not modify the leading behaviour, but in opposition 
to the Heisenberg model, it becomes relevant if 7(A^) = 7*(iV). In fact, the minimization 
of (^) gives 

L{N) ^ A^e^^^ . (30) 

Now, if 7(A^) = 'j*{N), the exponential dependence on A^ disappear^ and the polynomial 
dependence becomes relevant: therefore, a linear increase of the domain size is predicted 
(in the Ising model!) if Vex > N. 

We can summarize the main results in the following way. In both models (Ising and 
Heisenberg models), L has an exponential dependence on {J/2gw) ~ {Edw/g)- When the 
thickness is increased, in very rough terms g — > '~f*{N)g, but E^w has a different expression 



in the two models: Edw ~ J ^ l{N)J in the Ising one, and Edw ~ vJX —^ J'y{N)Edw in 
the Heisenberg one. So, in the limit 7 = 7*, the exponential dependence on A^ disapperas in 
the Ising model, but it remains in the other one. This dependence implies that L decreases 
with A^ in the Heisenberg model, whilst in the Ising model, if the exponential doesn't depend 
on A^, we have the pre-exponential dependence left: a dependence which makes L increase 
with A^ and which becomes a linear dependence in the continuum approximation. 

We can ask how the result for the Heisenberg model is modified, if the anisotropy is 
of surface-type: since the domain wall energy decreases with A^, narrower domains are 
expected! This result contradicts the findings on the closure domain structures [l^, based 



on the modified Kittel's theory. The reason is the following: in this theory, the anisotropy 
contributes to the total energy through the closure domains, where the magnetization lies 
in the film-plane, whilst the domain-wall energy is tacitly supposed not to depend on it. 

Instead, in our theory no closure domain is present —as expected, if the film is very 
thin— but the domain wall energy is proportional to the square root of the anisotropy. 

The conclusion of the present section is that —as a general rule— in a thin film domains 
are expected to get narrower upon increasing the thickness. This behavior has been found 



experimentally by Speckmann et al. [|T^ in ultrathin Co/Au(lll) films (above A^ = 4), by 



Bochi et al. |jT9[ in Ni/Cu/Si(001) films, and by Gehanno and Samson pO[ in ordered FePd 



thin films. So, a contradiction seems to exist between Ref. |T^ and Ref. |jT8|], both concerning 



Co/Au films. This contradiction might be solved guessing that the linear increase observed 



in Ref. WU accords with the weak increase of Ref. llTSi, oberved below A^ = 4. This behavior 



would correspond to the regime A^ < Vex, which gives —for an Ising model!— an increasing 



^In the limit iV ^ 1, the distinction between N and {N + 1) is not relevant. 
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L{N). Above four monolayers, we entry in the regime A^ > r^x (a narrow regime which is 
missed in Ref. ^^), which is characterized by a strong decreasing of L{N). 

Anyway, the previous picture is not totally satisfactory, because it assumes the validity 
of the Ising model, or —in other words— of domain walls of atomic thickness, whilst walls 
extend over hundreds of lattice constants. Furthermore, all the previous "magnetic" models 
assume perfectly flat and periodic arrays of spins, even if defects are known to play a major 
role in determing domain structures. For example, AUenspach and Stampanoni [17| report 



an increase of the average lattice constant, upon increasing A^. This means that dislocations 
exist in the film, and also that the exchange coupling constant is expected to decrease. Since 
J strongly determines the actual value of L, we can even suppose that the observed increase 
of -^(A^) is indeed due to the lowering of J, and therefore of the domain wall energy. 

4 Domain structures and reorientation transition 

In the part of the previous section which was dedicated to the Heisenberg model, we shew 
the possible different ground states of an ultrathin ferromagnetic film, subject to an easy- 
axis anisotropy and to the dipolar interaction. In Fig. ^ we show the phase diagram at 
T = 0, in the space of the parameters / = K/fi, L^ = J/Q. Upon decreasing the ratio 
K/Q, the system passes from a collinear ± state, to a collinear || state, across a domain 
structure. 

As we have pointed out, the _L state is never the ground state, in the thermodynamic 
limit. Anyway, from an experimental point of view, the question is different. First of all, 
domains may be even larger than the sample size, so that —in principle— the _L state may 
be the true ground state of a finite system. Furthermore, our estimates of the domain size 
were valid at T = 0, for a perfect sample. Whatever the domain size is, the _L state —which 
is metastable— is attained by the application of a magnetic field pulse: a typical procedure, 
before a magnetization measurement. Therefore, an important quantity is the relaxation 
time of the collinear state: if it is (much) longer than the measurement time, the _L state 
may be considered as the "ground state", whereas it will decay in a domain structure, in 
the opposite limit. 

Experimentally, the phase diagram in the parameter space may be "visited" by varying 
the thickness and/or the temperature. In fact, the easy-plane effect of the dipolar inter- 
action increases with A^, whilst the easy-axis anisotropy generally keeps constant, because 
of its surface character. This means that the "effective" ratio K/g decreases, and indeed 
a transition from the _L state to the || one is observed, when the thickness is augmented. 
This transition is called reorientation phase transition (RPT), and if a microscopic probe 
is used to measure the magnetization (for example, the spin polarized scanning electron 
microscopy), magnetic domains are observed in the "transition" region between the two 
collinear states |^ . 



The same phenomenon of reorientation may be observed upon increasing the tempera- 
ture, before the usual ferro-paramagnetic transition takes place. In this case, the theoretical 
study of the transition requires to take into account thermal fluctuations, which weaken the 
different coupling parameters {J,K,Q). The RPT is the result of a strongest decrease of 
the uniaxial anisotropy, with respect to the dipolar coupling: so, K{T) may be smaller 
than fi(T) —and therefore the || state the preferred one— even if K{T = 0) > Q{T = 0). 
The T-dependence may be obtained through a perturbation theory [^ ^, or through a 
renormalization group analysis [^]. We will not enter into details here, nor we will give the 
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Figure 3: "Phase diagram" at T = 0, in the plane {f^Ld). No units are shown on the / 
and Ld axis, because the diagram has only a qualitative meaning. The dotted line / = 1 
corresponds to the separation between the _L and || FM configurations. By taking into 
account domain structures, domains rise in all the region on the right of the full thin line, 
but they are observable only on the left of the dashed line. So, starting from a thin film 
which is perpendicularly magnetized, we may observe a reorientation phase transition upon 
increasing the temperature, or upon increasing the thickness. In a region around / = 1, 
domains of observable size appear. 



explicit dependences on T; however, we want to discuss an important feature of the RPT: 
the order of the transition. 

Firstly, let us address the question in the parameter space. If we consider the true 
ground state of the system, the answer is clear: the transition is continuous, because the 
domain structure (see Fig. la), in proximity to / = 1, evolves towards the || state, passing 
through a shrinking of the domain size and the continuous change of s (the maximal value 
of \Sz\ inside each domain) from s = 1 to s = 0. This picture clearly shows the importance 
of introducing such variational parameter in the study of domain structures. 

If domain structures are neglected, as in almost the whole body of theoretical papers on 
the RPT, the answer is more subtle, and —above all— it depends on details of the system. 
Limiting the investigation to coUinear configurations corresponds, in the _L phase, to study a 
metastable stateQ. Let us start with the simplest possibility: a single monolayer, subject to 
the Hamiltonian Ti = Hex+'Han+'Hdip- We have shown that Tidip has an effective easy-plane 
effect (TidJ^): as a consequence of this, the energy (per spin) of a collinear configuration 
writes: 

E{9) = -{K -^1)00^^9 , (31) 

where 6 is the angle between the magnetization and the z axis. It is manifest that E{6) allows 
only two possible minima, according to the sign of {K — Q): 6 = and 6 = 7i/2. So, in the 
parameter space, the transition is obviously discontinuous, because no canted configuration 



^It should be stressed that in the (very) narrow interval /„ 
locally stable. 



< / < 1, there is no collinear state which is 
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{0 < 9 < 7r/2) is allowed: so, it is not possible to change with continuity the value of 6 
which minimizes E{9), from = (_L state) to 6 = 7r/2 (|| state). Nevertheless, this result 
may be modified for two reasons: (i) the introduction of higher order anisotropies [^, and 
(ii) an increase of the thickness [p^j . 

In the former case, a term like K4Cos'^6 is added to E{6), and consequently E'{6) = 
allows solutions other than ^ = 0, tt/2. This means that the minimum is expected to 
change continuously from ^ = to ^ = tt/2, upon varying the parameters. In the latter case 
(A^ > 1), a possible canting of the magnetization is determined by the break of translational 
invariance in the perpendicular direction, due to the surface character of the easy-axis 
anisotropy |^6|], or even by the dipolar interaction itself, whose easy-plane effect varies a 
little from plane to plane flOl . 



From an experimental point of view, case (i) is much more important: in fact, the 
extension of the region where canting is really observed depends on K^/K in the first case, 
and on Q/J —or K/J— (two much smaller quantities) in the second one. In other words, 
variations of the effective anisotropy from layer to layer change the order of the RPT (from 
discontinuous to continuous), but the transition from _L to || keeps very sharp. 

Till now, we have analyzed the order of the RPT in the parameter space, at T = 0. 
What happens if the transition is induced by an increase in temperature? We will consider 
the simplest case: a single monolayer in absence of fourth-order anisotropies. In fact, this 
is the only case which shows a first order transition in the parameter space. 

We have said that spin-wave interactions lead to a renormalization of the coupling pa- 
rameters: (J, K, Q) -^ {J(T), K(T), Q{T)). Anyway, if the form of the Hamiltonian keeps 
the same, the energy of a coUinear configuration: E{6, T) = —{K{T) — ^{T)) cos^ 6 cannot 
have other minima than ^ = 0, it/2, and therefore no canted configuration is possible. This 



is the result of several authors [El, 12^, E^l |2S[, but it does not agree with recent calcula- 



tions by Chui [^ . This author claims that a single renormalization group equation for the 
dipolar coupling is not sufficient: under the flow of the renormalization group, the Hamil- 
tonian would acquire an additional term which changes the order of the transition, as the 
introduction of a quartic anisotropy made the transition a continuous one. Unfortunately, 



the lack of details in Ref. [29| prevents from drawing any conclusion. 



5 Domain structures at finite temperature 

We have seen how the competition between Hex and Tidip gives rise —at T = 0— to an 
ordered domain structure. This network of domain walls can be seen as a two dimensional 
solid, described by some positional and orientational order. This solid will be seen to 
support long wavelength excitations (wall meandering) as well as topological excitations 
(dislocations and stripe rotation domain walls). 

In the Ising model, a single domain wall can be considered the interface between regions 
of opposite magnetisation. In two dimensions, the free energy o"(T) of such an interface 
goes to zero at the critical temperature Tc of the Ising model; nevertheless, it represents 
a one dimensional interface, which is rough at any finite T. In the present case, "rough" 
means delocalized. This would imply that the solid is always in a fioating phase, with the 
walls which are depinned from the substrate. 

This is the scenario considered by Gehring and Keskin |3^, who simply replace the 
domain wall energy, by its free energy. Therefore L ^ exp [a{T)/g] decreases with T, till 
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L becomes of the same order of the domain wall size: from now on, no regular structure 
exists and a disordered phase sets in. 

Czech and Villain |^ had pointed out that a floating phase exists only at temperatures 
higher than Tr k. g. Indeed, neighbouring walls interact both entropically and via the long- 
range dipolar coupling. At sufficiently low T, entropic interactions can be neglected and 
a domain wall consists of straight pieces of wall separated by kinks. The typical distance 
between kinks is £ ~ e^"^/^, where (2 J) is the energy of a kink. A "coarse-grained" version of 
the system is obtained by sampling the m-th wall at the points yp = pi', where i' < i. The 
resulting dynamical variables are the displacements {«„} of the sampled points, with respect 
to the T = positions. The value i' is choosen sufficiently small so that \u^ — u^^\ < 1, 
but sufficiently large in order to suppose the different pieces of wall as independent. Even 
if kinks are absent {\u^ — u^^\ = 0), the rigid translation of a wall costs dipolar energy. 
The final expression for the free energy writes [^: 



m.p 



lll^m - < y + J±iul, - <-l) 



(32) 



which looks like an anisotropic Solid on Solid (SOS) model. The coupling constant in the 
wall direction (^|| = Tin (£/£')) is determined by the kink energy {= Tlni = 2J) and by the 
entropy S = Ini' of each piece of wall. On the other side, the coupling constant J± = gi'/L"^ 
is nothing but {d'^/dL'^){—glnL) times i'. 

It is easily seen that Jy ^ J± f^ g when T ^ g: this means that the roughening 
temperature of the solid of domain walls is not Tr = 0, but rather T^ k. g. At smaller 
temperatures, Czech and Villain suggest that L changes discontinuously (because of the 
discrete nature of the lattice), however keeping itself a too large value to be experimentally 
observable. At higher temperatures the network of domain walls represents a floating solid. 

We close this part by mentioning the mean-fleld result P|. In this approximation, the 
free energy writes: 



MF 



Y,(T-J{k))mkm.k + cTY.mt + ... , (33) 



where c is a numerical factor and J{k) is the Fourier transform of the exchange and dipolar 
coupling between spins. 

The value km which minimizes the quadratic form is of order (g/J), so that L{Tc) ~ J/g, 
whilst Tc = J (km) — J — 0{g). In other words, mean-fleld theory predicts relatively small 
domains at the transition (of size of order of the dipolar length) and that Tc is only slightly 
modifled by dipolar interactions (at least in the limit g <^ J). 

The most detailed and complete study on the different phases of a two dimensional 
network of domain walls has been done by Kashuba and Pokrovsky 0] and by Abanov, 
Kalatsky, Pokrovsky and Saslow ||31[ (hereafter called K2APS theory). They consider a 
Heisenberg model: because of the possibility of a Reorientation Phase Transition, the phase 
diagram is more complex. Anyway, if we neglect the phases where the order parameter is 
the magnetization, and we retain only the phase diagram of the two dimensional solid of 
domain walls, we obtain a system which does not differ too much from the Ising case. Indeed, 
in the limit L ^ w it is possible to replace the Heisenberg domains, separated by Bloch 
walls, with Ising domains, separated by empty walls. Spin-wave fluctuations are taken into 
account through a renormalization of the different coupling parameters. (This procedure is 
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expected to be correct if L is much larger than the length Lr at which renormalizations stop. 
Since this length is given by the relation J/ L\ k, A, the above condition equals L ^ w.) 

Once the spins have no in-plane component, the "volume" interaction in (1/r^) between 
domains can be reduced —through the Gauss's theorem— to an interaction in (1/r) between 
domain walls. This allows to pass from spin variables (Si) to "displacement" variables 
{u{x,y)) of the domain walls, and the energy of the system can be written as a functional 
of the field u{x,y). 

In very general terms, we can say that the energy must contain three terms: i) A 
compression term, expressing the response of the system to a changement of L. This term is 
proportional to {dxuY- ii) A bending term, expressing the energy due to a curvature {dyU) 
of the domain wall, iii) A stripe orientation term, expressing the existence of preferred 
directions for domain walls. This term, proportional to {dyuY, has an extreme importance 
because its presence drastically alters the behaviour of the system. 



By taking into account all the above mentioned terms, we obtain [|14|, ^ 
/ / dxdy 



E. 




(34) 



The first term is the "continuum" version of the term proportional to J±_ in the SOS Hamil- 
tonian by Czech and Villain [§]. In fact, according to K2APS theory, Ks ~ fi/L, which 
becomes, after the transformation x —>■ Lx, Ks ~ fl/L"^. This expression equals J±, once 
we have put i' = 1. 

We could therefore be induced to conclude that the third term corresponds to the one 
proportional to Jy in the same SOS Hamiltonian. Indeed, the question is more subtle, as 
suggested by the fact that z/ is not simply proportional to the domain wall energy. In fact, 
since 6 = dyU is the local orientation of the domain wall, in absence of anisotropic terms 
the energy Es can not explicitely depends on it. 

The microscopic nature of z/ is the following: once we have fixed all the parameters, 
the domain wall energy per unit length depends on the orientation of the wall itself. If 
we consider an Ising model on a square lattice, a domain wall in the direction (0,1) —or 
(1,0)— is less expensive that the same domain wall in the (1,1) direction: E^wiO, 1) = 2J/a 
and Edwi^A) = 2J/(a/v^) = \/2£'d^(0, 1), where we have exphcitly written the lattice 
constant a. So, on a square lattice we will have a stripe orientation energy with periodicity 
(7r/2). Its expansion in the vicinity of a minimum gives rise to the third term of Eg. 

The correct expression of u for the Heisenberg model is given in |14|: u = u^x = i; (^); 



where the factor (1/L) simply translates the continuum description of walls at distance L. 
We want to point out that v may have also a different origin: an in-plane anisotropy, for 
example. It is easily seen that in this case v = Van ~ K\\w/L, where K\\ is the single-ion 
in-plane anisotropy. So, Van > ^ex if (-^||/^) > {^1 J)i ^ condition which may be fulfilled in 
some systems. Furthermore, the periodicity of //„„ may be different from that of z/g^ (see 
the end of this section). 

Finally, the bending energy term (proportional to fi) —which was neglected in — 
is proportional to the square of the curvature (dyU), because of the required invariance 
u -^ —u. 

Once we have explained the origin of the energy Es-, let us discuss the T-dependent 
behaviour of our system. Two quantities are of interest: the mean square fluctuations of 
the position (u^), and the mean square fluctuations of the angular deviations {{dyuY). The 
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first ones are given by: 




KspI + vpl + ^lpl ^fK~s J yWT^ \ z/p2 + pip; 



K, 



(35) 



Since an ultraviolet cutoff is naturally supplied by the discrete nature of the system {\py\ < 
1), we are mainly concerned with a possible infrared divergence. When p^ — ;> 0, the integral 
has always a divergence: a logarithmic one if z/ 7^ and a power-like one if z/ = 0. The 
difference is relevant, because the calculation of the correlation function {u{Q)u{r)) shows a 
power-like decay, and therefore the existence of a quasi long-range order (QLRO), solely in 
the first case: in absence of a stripe orientation energy, domains are completely delocalised 
and the derivative {dxu) has even no meaning. The calculation of {{dyuY) simply needs 
to add a factor {p^) at the numerator of the previous integral: because of this factor, the 
angular deviations are always finite, even if z/ = 0. 

By summarizing: in the "real" case z/ 7^ 0, there is an orientational long-range order 
(LRO) and a positional QLRO. A magnetic domain phase with these properties has been 
indeed observed by AUenspach and Bischof [^] in an ultrathin film of Fe/Cu(100), with 



stripe domains oriented along a high symmetry direction in the plane. 

The absence of LRO for the position of domain walls apparently clashes with the results 
by Czech and Villain. Anyway, no real contradiction exists: at the extremely low temper- 
atures considered in |]^ (T < fi), excitations in the Ising model are given by kinks, and a 
discrete description of domain walls is necessary. In this case, we obtain the SOS Hamil- 



tonian (|32D , to which corresponds a finite depinning temperature. Conversely, the picture 
given by K2APS theory is a continuous one, valid at not too low temperatures (T > fi), 
where even the Czech and Villain's theory foresees a floating phase. 

Let us come back to Eq. (^4]): at sufficiently large lengthscales, the higher order term 
{dyU)"^ is negligible, and Es keeps an anisotropic quadratic form. This energy is well 



known |5^ to correspond to an algebraical decay of the fluctuations of the order parame- 
ter. The corresponding solid is called a smectic crystal by K2APS: the QLRO is destroyed 
by the proliferation of pairs of dislocations|] and the corresponding critical temperature 
writes according to the Kosterlitz-Thouless expression |3^ (P stands for "positional"): 
Tp = ^K'lP' = ^ ^^'^ L'^, where K' = \/Ksv is the (renormalized) effective elastic constant 
and 6 = 2L is the Burger vector corresponding to the dislocation. 

Above Tp the compression coefficient vanishes and domain walls are delocalized. What 
happens to the orientational order parameter (hereafter explicitely written through the angle 
6')? The bending energy rewrites as ^{dyiOY ^ whilst for the orientational stripe energy 
the complete expression will be used: y|(1 — cos46'). We must also add a term of the 
form ^{dx'OY^ which would correspond to f (c^^j^m)^ in the "old" coordinates and which was 



V 



' K fdeV fi fde\' . .^ 



(36) 



irrelevant to determine the large scale properties of Es- Finally, we will have [^ 
E^ = dx'dy 

where {x',y') are local coordinates. 

^A pair of dislocations of opposite Burger vectors rises from the "collapse" of a piece of two neighbouring 
walls. 
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An important point is the sign of n, which may be positive or negative: its mean-field 
calculation gives Hmf ~ —^L, but both thermal fiuctuations and corrections due to free 
dislocations give [^ a positive contribution to k. An orientational ordered phase will exist 



at T > Tp, only if k > 0. In this case, because of the anisotropic term proportional to u 
—which assures an orientational LRO— the order parameter 6 has a discrete Z2 symmetry: 
domain walls, in a square lattice, can be directed either along a given direction, or along the 
direction perpendicular to it. For this reason, this phase has been called Ising nematic phase. 
At temperatures slightly larger than Tp, the two directions are not equally probable and 
the Z2 symmetry is broken. Two contiguous regions where domains are directed along the 
two different "easy-axes" , are separated by a so-called stripe rotation domain wall, studied 
in detail by Abanov et al. |^. At a temperature To > Tp we therefore expect an Ising- 



type transition towards a liquid phase. Since there will be two (mutually perpendicular) 
directions which are more likely, this phase is called a "tetragonal liquid". In this phase, 
square domains should be still observable, with a magnetization inside each domain which 
is perpendicular to the surface. This phase may evolve towards a planar phase, via a 
Reorientation Phase Transition, or directly towards a paramagnetic phase. 

For a negative k, the Ising nematic phase does not exist and there is |^ a first-order 
phase transition between the (smectic) crystal and the (tetragonal) liquid. 

Till now, we have supposed that the anisotropic stripe energy has a (7r/2)-periodicity 
and therefore domains orient preferentially along two perpendicular directions. We can 
suppose that in certain cases the in-plane symmetry is reduced, for example because of an 
uniaxial distorsion, which determines a quadratic spin anisotropy (which looks like Ti'l^ = 
— -^ J dxSy] or a modification of Hex such that the exchange coupling J differs in the two 
spatial directions x and y. In both cases, the existence and the nature of the smectic- 
nematic transition at T = Tp are not changed, because for T < Tp only the expansion at 
small 6 of the stripe energy is relevant, not its global periodicity. Conversely, for T > Tp 
we are interested in the orientational order, and for a vr-periodicity the resulting scenario is 
completely different. In fact, in this case there is only one preferred direction for domain 
walls and no symmetry breaking transition may take place at T = To. The system behaves 
as a nematic in an electric field. The orientational order is not destroyed, till a "magnetic" 
transition towards a planar phase, or a paramagnetic phase, occurs. 

It is interesting to compare the given picture with the results of a Monte Carlo simulation 
on an Ising square model, due to Booth et al. [3l|. The authors study a 64 x 64 spins sample 



with a ratio J/fi ~ 9, such that L ~ 8 at T = 0. They single out three different phases, 
through the existence of two maxima in the specific heat C. By the comparison with 
Monte Carlo configurations at different temperatures, they conclude that there is (I) A 
low-T phase with an orientational order of domain walls; (II) A middle-T phase with no 
orientational order, but still well defined domains; (III) A high-T phase, which is fully 
disordered. 

Authors suggest that (I) resembles the smectic crystal phase, and (II) the tetragonal 
liquid one. Since their model can not undergo a RPT, (III) is the paramagnetic phase. 
Authors also add that such parallel suffers a weak point: the transition between (I) and (II) 
seems to be continuous, whilst the transition between the smectic crystal and the tetragonal 
liquid would be first-order, or there would be a further phase (the Ising nematic one) in 
between. 

In this respect, some comments are in order. Authors |Q study only the orientational 



order, and no details are given about the positional one. Besides, domain size is rather 
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small and dislocations are not visible in the figures reported in ||3^. However, the most 
important point is the nature of the transition between the smectic crystal and the Ising 
nematic phase. This transition corresponds to the loss of the positional QLRO and has 
the nature of a Berezinsky-Kosterlitz-Thouless transition. Such transition, linked to the 
unbounding of pairs of dislocations, is not "visible" through the specific heat C, because it 
corresponds to an unobservable essential singularity and C has no maximum for T = Tp. 

We argue that the observed phase transition between (1) and (11) should correspond to 
the transition between the Ising nematic phase and the tetragonal liquid. 

6 Conclusions 

In this paper we have reviewed theoretical results on magnetic domains in ultrathin films. 
Two-dimensional domain structures have a more fundamental meaning than the correspond- 
ing ones in bulk systems: in fact, in three dimensions domains strongly depend on the shape 
of the sample, and in a given sample are generally present for any direction of the magne- 
tizationQ. Conversely, in two dimensions shape effects are negligible in the thermodynamic 
limit, but domains are absent if the magnetization lies in the film plane. Furthermore, the 
interplay between uniaxial anisotropy and shape anisotropy possibly determines a reorien- 
tation transition: domains are expected to appear just before such transition takes place 
(upon increasing the thickness and/or the temperature). 

The main limit of the present article is to have neglected the possible application of a 
magnetic field. Because of the "reduced" quantity of matter, in a thin film the energy gain 
of a domain structure with respect to a collinear state is extremely low (see Fig. lb). As a 
consequence of this, a lower magnetic field is necessary to obtain a single domain state. 

A theoretical comprehension of the T-dependence of the network of domain walls has 
received an important impulsion with the K2APS theory, but the possibility to discover 
in a real system the proposed sequence of phase transitions is limited by the smallness 
of the region of temperatures where domains are observable. Furthermore, comparison of 
theoretical predictions with experimental data is made more difficult in ultrathin films by 
the increased importance of defects (dislocations and vacancies) and by the roughness of 
the interfaces. 

Finally, simulation data start to be available, but they often concern Ising systems, and 
domain sizes are generally far from being realistic. 
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